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Abstract 



In this paper, we provide a detailed description of our recent analysis and 
determination of the frame-dragging effect obtained using the nodes of the 
satellites LAGEOS and LAGEOS 2, in reply to the paper "On the reliability 
?h ! of the so-far performed tests for measuring the LenseThirring effect with the 

LAGEOS satellites" by L. Iorio (doi: 10.1016/j.newast.2005.01.001). First, we 
discuss the impact of the Jm uncertainties on our measurement and we show 
^ ■ that the corresponding error is of the order of 1 % of frame-dragging only. 

We report the result of the orbital simulations and analyses obtained with 
and without J4 and a with J4 equal to its EIGEN-GRACE02S value plus 12 
times its published error, i.e., a J 4 equal to about 611 % of the value adopted 
in EIGEN-GRACE02S, that is j 4 = 6.11 x (-1.41 • lO" 11 ) ~ -8.61 ■ 10" 11 . 
In all these three cases, by also fitting the final combined residuals with a 
quadratic, we obtain the same value of the measured Lense-Thirring effect. 
This value differs by only 1 % with respect to our recent measurement of 
the Lense-Thirring effect. Therefore, the error due to the uncertainties in J2 n 
in our measurement of the gravitomagnetic effect can at most reach 1 %, in 
complete agreement with our previously published error budget. Our total 
error budget in the measurement of frame-dragging is about 5% of the Lense- 
Thirring effect, alternatively even by simply considering the published errors 
in the J2 n and their recent determinations we get a total error budget of the 
order of 10 %, in complete agreement with our previously published error 
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budget. Furthermore, we explicitly give the results and plot of a simulation 
clearly showing that the claim of Iorio's paper that the J2 n uncertainty may 
contribute to up a 45 % error error in our measurement is clearly unsubstanti- 
ated. We then present a rigorous proof that any "imprint" or "memory" effect 
of the Lense-Thirring effect is completely negligible on the even zonal harmon- 
ics produced using the GRACE satellites only and used on the orbits of the 
LAGEOS satellites to measure the frame-dragging effect. In this paper we do 
not discuss the problem of the correlation of the Earth's even zonal harmonics 
since it only refers to our previous, 1998, analysis with EGM96 and it will 
be the subject of a different paper; nevertheless, we stress that in the present 
analysis with EIGENGRACE02S the total error due to the static Earth grav- 
ity field has been calculated by pessimistically summing up the absolute values 
of the errors due to each Earth's even zonal harmonic uncertainty, i.e., we have 
not used any covariance matrix to calculate the total error but we have just 
considered the worst possible contribution of each even zonal harmonic uncer- 
tainty to the total error budget. We also present and explain our past work 
on the technique of measuring the Lense-Thirring effect using the LAGEOS 
nodes and give its main references. Finally we discuss some other minor points 
and misunderstandings of the paper by Iorio, including some obvious mistakes 
contained both in this paper and in some other previous papers of Iorio. In 
conclusion, the criticisms in Iorio's paper are completely unfounded and mis- 
directed: the uncertainties arising from the possible variations of J2„ are fully 
accounted for in the error budget that we have published. 



1 Error due to the J 2n in the 2004 measurement of the Lense- 
Thirring effect 

In order to discuss the error analysis and the total error budget of our measure- 
ment of the Lense-Thirring effect (Ciufolini and Pavlis 2004), we first stress 
that in the data reduction of our recent measurement of the Lense-Thirring 
effect (Ciufolini and Pavlis 2004) we have used the value of J 4 = —1.41 • 10~ n , 
adopted by GFZ in the EIGEN-GRACE02S Earth gravity model (Reigber et 
al. 2005), and we have fitted our combined residuals with a secular trend only 
plus a number of periodical terms. We can of course introduce J4 as a free pa- 
rameter in our fit (see below). In this case, together with the measurement of 
the Lense-Thirring effect, we also measure the effect of the secular variations 
of J 2 , J4 and J 6 on the combination of the nodal longitudes of the LAGEOS 
satellites; this is described by a jf^ ective (Eanes 1995) in our combination, 
which includes the effect of the secular variation of the higher even zonal har- 
monics. In (Ciufolini et al. 2005) it is indeed reported an effective value of 
jEffective ^ _i 5 . io~ n for the combination of the LAGEOS satellites nodes, 
which is consistent with the EIGEN-GRACE02S model since, on our combina- 
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Fig. 1. Fit of the residuals of the nodes of LAGEOS and LAGEOS II, using our 
combination and the Earth's model EIGEN-GRACE02S, with a secular trend plus 
six periodical terms. The slope is fi ~ 0.99 (Ciufolini and Pavlis 2004). 

tion of the nodal longitudes of the LAGEOS satellites, it just represents a 6 % 
variation of the value given with EIGEN-GRACE02S and, however, it includes 
the effect of any higher J2 n , with 2n > 4; this value is also fully consistent 
with our published result of a Lense-Thirring drag equal to 99 % of the general 
relativity prediction with uncertainty of 5 % to 10 %. It is easily seen, even by 
visual inspection, that our combined residuals would clearly display any such 
quadratic term. Indeed, in figure (1) we show the residuals obtained using 
the J4 value given with EIGEN-GRACE02S that should be compared with a 
simulation of the orbital residuals, shown in figure (5), obtained using in the 
data reduction a strongly unrealistic value of J 4 corresponding to the value 
adopted in EIGEN-GRACE02S plus 12 times its published error and which 
produces a ~ 45% variation of the secular trend as claimed by the author of 
(Iorio 2005)! It is clear that only the first figure can be simply described by a 
linear dependence. 

In the EIGEN-GRACE02S model (Reigber et al. 2005), obtained by the 
GRACE mission only, the Earth gravity field was measured during the period 
2002-2003. Corrections due to J 2 and J 4 were then applied to this 2002-2003 
measurement in order to obtain a gravity field model antecedent to 2002- 
2003. These values of J2 and J4, used by the GFZ team, are J2 = —2.6 ■ 10 -11 
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and J4 = —1.41 • 10 11 and they were measured on the basis of completely 
independent 30-year observations before 2002. 

Let us describe the result of the orbital analyses using the orbital estimator 
GEODYN with and without a contribution of J4 = —1.41 • 10~ n ; and the 
result of a simulation with a J4 equal to 611 % of the value given in EIGEN- 
GRACE02S, i.e. J 4 ~ -8.61 • 10" 11 . 

First, we stress again that in the case of not applying the J 4 EIGEN-GRACE02S 
correction to the orbital analysis, its effect can be clearly identified (but nev- 
ertheless completely fitted for, see below) even by visual inspection since it 
shows up as a hump in the combined residuals, whereas in our fit of (Ciufolini 
and Pavlis 2004) shown in Fig. 1 the absence of any such quadratic effect is 
obvious. Then, since the effect of the time variation J 4 shows up as a quadratic 
effect in the cumulative nodal longitude of the LAGEOS satellites, the com- 
bined residuals of LAGEOS and LAGEOS 2 may be fitted with a quadratic 
curve, together with a straight line and with the main periodic terms. Thus, 
by fitting the raw residuals obtained without any J 4 in a test — analysis, see 
Fig. 2, we measured a j^ff ectlve ^ — 1.5 • 10~ n , which includes the effect of j 6 
and of higher even zonal harmonics on the combination of the LAGEOS nodes. 
On other hand, by fitting the combined residuals obtained with the EIGEN- 
GRACE02S correction of J4 = — 1.41 • 10~ n , we measured a §j^f ecttve Q f less 
than — 0.1 ■ 10~ n , see Fig. 3, and this fitted difference is in complete agreement 
with the previous case. In all the cases of inclusion of an anomalous J4 in the 
data analysis, we obtained a measured value of /i = 0.98, see Figures 3 and 4. 

Therefore, the small value of the unmodelled quadratic effects in our nodal 
combination due to the unmodelled J 2n effects (with 2n > 4) corresponds to a 
change in the measured value of the Lense-Thirring effect of about 1 % only. 
As a consistency test, using the value j^ff ecttve — — 1.5 . 10 -11 , obtained from 
fitting the combined residuals (which is only about 6 % larger than the value 
j 4 = -1.41 • 10" 11 given in the EIGEN-GRACE02S model), we have again 
generated the residuals and fit these with a straight line only plus the main 
periodic terms. It turned out that the change of the measured value of frame- 
dragging was about 1 % only with respect to the case of using J 4 — 1.41 • 10~ n 
and fitting the residuals with a straight line plus the main periodic terms 
(Ciufolini and Pavlis 2004), i.e., this test resulted in a measured value of fi = 
0.98. The result of this test agrees with the case of fitting the combined orbital 
residuals with a straight line plus a quadratic curve, indeed in all cases we got 
= 0.98, see Figures 2, 3 and 4. As a third case, using in the data reduction a 
highly unrealistic value of J 4 equal to the value adopted in EIGEN-GRACE02S 
plus 12 times its published error, i.e., a J4 equal to about 611 % of the value 
given in EIGEN-GRACE02S, that is j 4 = 6.11x(-1.41-10 _11 ) ~ -8.6M0 -11 , 
by fitting the residuals with a straight line plus the main periodic terms plus 
a parabola, we again obtained the same measured value of frame-dragging, 
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Fig. 2. Fit of the residuals of the nodes of LAGEOS and LAGEOS II, using our 
combination and the Earth's model EIGEN-GRACE02S, with a secular trend plus 
six periodical terms plus a quadratic. The residuals have been obtained without in- 
clusion of any J4 value in the orbital analysis, i.e., using J4 = 0. The Lense-Thirring 
slope is fj, ~ 0.98. 

i.e., \x = 0.98, see Fig. 4. Therefore, also in this case the measured value of 
frame-dragging differs by about 1 % only with respect to the case of using 
J4 = —1.41 ■ 10~ n without fitting a quadratic curve, i.e., the result reported 
in (Ciufolini and Pavlis 2004) and in Fig. (1). 

These analyses clearly show that even in the cases of - 100 % and + 511 % 
variations between the value of J4 of Mother Nature and the value used in 
the orbital estimation of (Ciufolini and Pavlis 2004), it is possible to fit for 
the jEff ectlve e ff ec t anc [ g e t back from the fit most of the simulated variation 
in J4 but especially it shows that the measured value of the Lense-Thirring 
effect can only be affected at the 1 % level by the J2 n uncertainties. Indeed, 
by fitting with a parabola the combined orbital residuals in these three cases 
corresponding to different values of J4 (zero, the EIGEN-GRACE02S value, 
and 611 % of this EIGEN-GRACE02S value), we always obtained \i = 0.98 
(Figures 2, 3 and 4, whereas in (Ciufolini and Pavlis 2004) we obtained, 
without fitting any parabola, fi = 0.99. 

In conclusion this 1 % variation (obtained by fitting the residuals including 
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Fig. 3. Fit of the residuals of the nodes of LAGEOS and LAGEOS II, using our 
combination and the Earth's model EIGEN-GRACE02S, with a secular trend plus 
six periodical terms plus a quadratic. The residuals have been obtained using the J4 
value given in the EIGEN-GRACE02S model. The Lense-Thirring slope is \i — 0.98. 

a parabola) with respect to the value of \x = 0.99 measured in (Ciufolini and 
Pavlis 2004) (without any parabola in the fit) gives the estimated error due 
to the J2n uncertainties in our measurement of /1, as reported in (Ciufolini 
and Pavlis 2004). In (Iorio 2005), it is claimed that the 5 % to 10 % error 
budget of our recent measurement has been strongly underestimated because 
of the J4 and Jq errors; indeed, in this paper it is claimed a possible error 
as large as 45 % due to the J 2n uncertainties. However, this statement is 
simply nonsense based on what we explained in the previous paragraph and 
only shows a lack of understanding of both the real effect of the J2 n , with 
In > 4, on the combination of the nodes of LAGEOS and LAGEOS 2 and of 
the technique of the least-squares fit. 

As already pointed out, it is critical to stress that anyone can immediately 
rule out a large error of this type even by visual inspection of Fig. 1, which 
is the fit of the combined orbital residuals of (Ciufolini and Pavlis 2004). A 
change, for example, of 45 % in the fit due to the J 2n effects (as claimed in 
(Iorio 2005)) would correspond to a superposition of a gigantic hump, with 
a height of several hundreds of milliarcsec, to our raw residuals that instead 
clearly show a straight line trend with a rate of 47.9 milliarcsec (compare Fig. 
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Fig. 4. Fit of the residuals of the nodes of LAGEOS and LAGEOS II, using our 
combination and the Earth's model EIGEN-GRACE02S, with a secular trend plus 
six periodical terms plus a quadratic. The residuals have been obtained using a J4 
value equal to 611 % of the j 4 value given in the EIGEN-GRACE02S model. The 
Lense-Thirring slope is \i ~ 0.98. 

1 and Fig 5)! 

Indeed, in Fig. 5 we have shown the result of a simulation using, as an example, 
a value of J4 = —8.61 • 1CT 11 , corresponding to the value given in the EIGEN- 
GRACE02S model plus 12 times the error of J4 given in (Cox and Chao 2002) 
(see below), for example with x-intercept in the year 2000. Not only, in this 
case, to get a 45 % error we would need an error in J4 that is 1200 % of what 
estimated in ref. (Cox and Chao 2002), but any such error would produce a 
huge deformation of the final residuals that can be easily identified even by 
visual inspection of the combined residuals and, however, completely fitted 
for using a quadratic as shown in Fig. 4 and as previously explained. A 45 
% error in the measured value of \i due to a superimposed quadratic effect 
looks as a clear nonsense! In other words, by comparing Fig. (5) with Fig. (1) 
(corresponding to our measurement (Ciufolini and Pavlis 2004)), it is clear that 
what is missing and misunderstood in (Iorio 2005) is that the only relevant 
secular effect that can mimic the Lense-Thirring effect is due to the errors 
in the values of the static J 2n coefficients used in our analysis, whereas the 
Jm will show up as quadratic effects and therefore large values of Jm can be 
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Fig. 5. Simulation of our combined residuals of the nodes of LAGEOS and LAGEOS 
II using our combination. The residuals have been obtained using a highly unrealistic 
value of j 4 equal to 611 % of the j A value given in the EIGEN-GRACE02S model 
with x-intercept in 2000. The Lense-Thirring slope, fitting only for a secular trend, 
is fi ~ 0.55, i.e., it has a change of about 45 % with respect to the fit of Figure 2. The 
quadratic effect in the combined residuals due to J4 can be clearly identified even by 
visual inspection. This plot should be compared with Fig. 1 (the result of (Ciufolini 
and Pavlis 2004)) and with Fig. 4 where these simulated combined residuals have 
been fitted by including a quadratic term in order to fit for the j^Sf ectlve effect. 

clearly identified over our long period of 11 years and can be fitted for, using 
a quadratic curve. For example, in the single, i.e., not combined, residuals of 
the node of LAGEOS and of the node of LAGEOS 2 of Fig. 6 and 7 is clearly 
observed a long term variation in the trend that, since it disappears in the 
nodes combination of Fig. 1, can be identified as an anomalous increase in the 
Earth quadrupole moment, corresponding to the effect observed by Cox and 
Chao (Cox and Chao 2002). 

Finally, since the measurements of J 2n and J 2n are completely independent 
(as previously explained), the corresponding errors are independent and this 
allowed us to take the RSS of these errors. Nevertheless, even by just adding 
up these errors we would get an error budget of 6.6 %, well within our quoted 
range of 5 % to 10 %. 
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. 6. Residuals of the nodal longitude of LAGEOS using EIGEN-GRACE02S . 




7. Residuals of the nodal longitude of LAGEOS 2 using EIGEN-GRACE02S. 
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Another way to look at the interconnection between the Lense-Thirring effect 
and the value of J4 is through a \ 2 plot, as shown in Fig. 8. The origin of 
this graph corresponds to the best fit of our combined residuals by including 
a parabola in the fit (i.e., it corresponds to j^f ective ^ -1.41 • KT 11 - 0.1 • 
10~ n = —1.51 • 10~ n and Lense-Thirring effect = 98% of the general relativity 
prediction) and the ellipse represents the confidence level (Press et al. 1986) 
for a 99 % probability. It has to be stressed that this ellipse of confidence is 
tentative in the sense that we have assumed an average tentative standard 
deviation for each combined residual; the full analysis using the covariance 
matrix obtained by the GEODYN data reduction is presented in (Ciufolini 
et al. 2005) . In this figure we have also included the value of J4 adopted 
in the EIGEN-GRACE02S model (Cox and Chao 2002). The two parameters 
are indeed correlated, but this figure shows that the maximum uncertainty in 
the Lense-Thirring parameter arising from the possible variation of J4 is just 
~ ±0.75% of its value. The complete error analysis published in (Ciufolini 
and Pavlis 2004) and (Ciufolini et al. 2005) quotes a larger total uncertainty 
between 5 % and 10 % since it also takes into account all the other error 
sources. 

Let us estimate with a different method the maximum conceivable error due to 
the the J2 n in our measurement of the Lense-Thirring effect. First we stress that 
our error budget relies on the validity of the EIGEN-GRACE02S model and 
of the GRACE mission. Questioning our error budget is basically equivalent 
to question the validity of the GRACE mission itself. 

As already pointed out, the Earth model EIGEN-GRACE02S was produced 
by GRACE measurements only, during the period 2002-2003, whereas the 
independent evaluation of J4 was produced on the basis of completely inde- 
pendent 30- year observations before 2002. The only constraint (based on the 
validity of the GRACE measurements) is that the J 4 correction applied to J 4 , 
must of course produce in 2002-2003 the same value of J4 that was measured 
by GRACE in 2002-2003, at least within the EIGEN-GRACE02S uncertain- 
ties in J4. Running then the orbital estimator GEODYN with or without a 
contribution corresponding to the estimated error in J 4 of ±0.6 • 10~ n (Cox 
and Chao 2002) and by fitting the residuals without a parabola resulted in a 
change of the corresponding measured value of frame-dragging of ± 7.29 %. 
This was verified in two ways: using GEODYN and an independent program 
written specifically to evaluate the J4 effect. By changing the observational pe- 
riod the measured value of frame- dragging would only slightly change (see the 
slope of Fig.l over subsets of the total observational period of 11 years) and 
the corresponding error would consistently change too, increasing over peri- 
ods farther away from 2002-2003. Nevertheless, over our observational period, 
the maximum J4 error was 7.29 %. Similarly, by fitting the residuals without 
a parabola, an error of 0.5 • 10~ n in J 6 would produce a 2.98 % error over 
our observational period. Then, by adding these errors to the error due to 
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Fig. 8. Confidence Ellipse. On the x axis is the Lense-Thirring effect; the origin, at 
x = 1, corresponds to the fitted value of frame-dragging. On the y axis is jEff ecUve .^ 
the origin, at y = 1, corresponds to the fitted value j^ff ecttve ~ —1.51 • 10 -11 . 
This confidence ellipse is tentative since we have assumed an average tentative 
standard deviation for each combined residual. This figure shows that the maximum 
uncertainty in the measured Lense-Thirring parameter arising from the possible 
variation of J4 is just ~ 0.75% 

tides we get a maximum error of 11 % and finally considering all the other 
independent error sources (see the Supplementary Discussion of ref. (Ciufolini 
and Pavlis 2004)) we get a total uncertainty of 12 % in our measurement of 
the Lense-Thirring effect, in substantial agreement with the maximum error 
of about 10 % reported in (Ciufolini and Pavlis 2004). Nevertheless, it is crit- 
ical to perform this worst case error analysis (without fitting our residuals 
with a quadratic) by using the most recent 2004 measurements of J4 and J§ 
(Cheng and Tapley 2004), see also section 2. Recently, a solution for the time 
variation of the even zonal harmonics up to degree 6 was determined using 
Satellite Laser Ranging data from 8 geodetic satellites over a time period of 
28 years from 1976 to 2003. The secular zonal variations for the degree up to 
6 have been determined using the existing long term Satellite Laser Ranging 
data from multiple satellites, including Starlette, Ajisai, Stella, Lageos I and 
II, Etalon I and II, and BEC using both long and short arc orbital analysis 
technique. The results of (Cheng and Tapley 2004) are J4 = —1.99 • 10~ n and 
Jq = 0.02 ■ 10~ n . Then, we have used these values of J4 and Jq and the or- 
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bital estimator GEODYN, according to the worst case analysis just described 
(without any fit for a quadratic term), to evaluate the impact of these differ- 
ent values of J 4 and J 6 in our measurement of the Lense-Thirring effect. We 
have then found a total difference of about 6.9 % in our measurement of the 
Lense-Thirring effect due to these new values of J 4 and J 6 (a precise statis- 
tical analysis based on the scatter of the various determinations of the J2 n is 
presented in a detailed paper in preparation (Ciufolini et al. 2005)). Then, by 
adding this error to the error due to tides we get a maximum error of less than 
7.9 % from time variations of the Earth gravity field and finally considering 
all the other independent error sources we get a total uncertainty of 9.3 % 
in our measurement of the Lense-Thirring effect, in full agreement with the 
maximum error of about 10 % reported in (Ciufolini and Pavlis 2004)- 

However, as explained above, we stress that the effect of such an error in 

jEffective ^ ^ obgerved 

in our combined residuals, nor is of course observed 
a larger error of this type (Iorio 2005) that would in fact produce in our 
combined residuals a huge parabola superimposed to the frame- dragging trend 
(see, for example, Fig. 5). Indeed, as previously explained, by fitting our raw 
residuals including a quadratic curve we have that the total effect of the J 2ri , 
with 2n > 4, on the combined node residuals corresponds to a total error in 
our determination of frame-dragging of about 1 % only. 

It is finally worth to stress once again that the measured value of the Lense- 
Thirring effect, /i = 0.98, obtained by fitting with a straight line and six 
periodic terms plus a parabola the raw combined residuals of (Ciufolini and 
Pavlis 2004), fit shown in Fig. 2, fully agrees with the measured value of 
the Lense-Thirring effect, fi = 0.98, obtained by fitting with a straight line 
and six periodic terms plus a parabola the raw combined residuals simulated 
over our observational period of 11 years using GEODYN without inclusion 
of any J 4 , i.e. with J 4 = (Fig. 3), and with a highly unrealistic value of 

j 4 = 611 % jflGEN-GRACEmS ^ 4) _ 

In conclusion, the 5 % uncertainty in our measurement of the Lense-Thirring 
effect is fully confirmed in spite of the serious misunderstandings of our anal- 
ysis in (Iorio 2005). Indeed, a number of further critical remarks are reported 
in (Iorio 2005) that show that his author did not carefully and seriously read 
our papers but only superficially read them and deeply misunderstood them. 

In (Iorio 2004) is written "The value j 4 ^ = —1.5 - lO^ 11 year" 1 measured with 
the combination of Eq. (5) is affected not only by J 6 and the other higher 
degree even zonal harmonics but, more importantly, by the LenseThirring sig- 
nature itself. Indeed, the combination of Eq. (5) is designed in order to only 
disentangle Jq and the LenseThirring effect. So, it is not admissible to use 
the so obtained j 4 , which is coupled by construction to the LenseThirring 
effect, in order to reliably and correctly measure the LenseThirring effect it- 
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self" This statement is completely false because the Lense-Thirring effect was 
measured in (Ciufolini and Pavlis 2004), as very clearly written in (Ciufolini 
and Pavlis 2004), by only fitting with a straight line (plus the main periodical 
terms) the residuals produced with the standard value of J 4 adopted in the 
EIGEN-GRACE02S model: no effective J 4 and no quadratic fit were used in 
order to obtain our result in (Ciufolini and Pavlis 2004). The quadratic fit, 
together with a straight line and harmonics terms, was only used as a "tool", 
to evaluate the strength of any unmodelled J 4 effect (as above explained in 
details). Furthermore, as clearly explained above, the statement "the so ob- 
tained j e J* , which is coupled by construction to the LenseThirring effect" is 
completely void of any possible meaning. 

Furthermore, this paper (Iorio 2004) is not properly referring to all the funda- 
mental papers published since 1984 on the use of the nodes of the LAGEOS 
satellites in order to measure the Lense-Thirring effect (see section 4). 

In the conclusion of (Iorio 2004) is written: "For a given observational time 
span and a given background reference Earth gravity model, two time series 
of the combined residuals, built up with and without the aforementioned J; 
GTR-free values in the background reference models, it should be analyzed. 
The difference between the so obtained fi parameters can, then, be evaluated." 
Again, in an official rebuttal letter of Ciufolini and Pavlis, dated November 
2004 (Ciufolini and Pavlis 2004), to answer one of the papers of criticism sent 
by Iorio for publication (rejected by the referees), the method to evaluate 
the impact of the J 2n uncertainties was clearly written and explained: " ... 
Running then the orbital estimator GEODYN with or without a contribution 
corresponding to the estimated error in J 4 of 0.6 • 10~ n (Cox and Chao 2002) 
resulted in a change of the corresponding measured value of frame-dragging 
by 7.29 %. This is exactly one of the methods of error analysis reported 
here. 

However, the most profound misunderstanding in (Iorio 2004) is the following 
"... If one is interested in extracting from such a signal a linear trend, it is not 
possible also to fit a parabolic noise signal without corrupting or distorting 
in some way the genuine linear trend of interest. Exactly in the same way, 
it is not possible to fit and remove the noise linear signal induced by J 4 , J 6 , 
. . . It is only possible to evaluate, as more accurately and realistically as 
possible, their impact on the measurement of the linear trend.". In a previous 
similar paper of the same author (Iorio 2004) is written: "The time series of 
the combined residuals of the nodes of LAGEOS and LAGEOS II consists of 
the Lense-Thirring linear trend as predicted by GTR, which has, in fact, to be 
measured, the linear trend induced by the static part of the uncancelled even 
zonal harmonics of the geopotential J 4 , J$ the parabolic signal induced 
by J 4 , Jq and the time-varying harmonic perturbations of gravitational and 
non gravitational origin. If one is interested in extracting from such a signal 
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a linear trend, modelled as ao + ait, it is not possible also to fit a parabolic 
noise signal, modelled as bo + b±t + b 2 t 2 , without corrupting or distorting in 
some way the genuine linear trend of interest. 

However, in the above mentioned official letter of Ciufolini and Pavlis to Io- 
rio, dated November 2004 (Ciufolini and Pavlis 2004), it was clearly writ- 
ten: "First, our error budget relies on the validity of the EIGEN-GRACE02S 
model and of the GRACE mission. Questioning our error budget is equivalent 
to questioning the validity of the GRACE mission itself. We stress that the 
Earth model EIGEN-GRACE02S was produced by GRACE measurements 
ONLY, during the period 2002-2003, whereas the INDEPENDENT evalua- 
tion of J4-dot was produced on the basis of completely independent 30-year 
observations BEFORE 2002. The only constraint (based on the validity of the 
GRACE measurements) is that the J 4 -dot correction applied to J 4 , MUST of 
course produce in 2002-2003 the same value of J 4 as that was measured by 
GRACE in 2002-2003." 

This means that one is not free to fit for a trend b\t due to J 4 in addition 
to the trend a\t but, on the basis of the 2002-2003 GRACE measurement 
of J 4 , any change J 4 , when applied to the EIGEN-GRACE02S static model, 
must of course be such to agree in 2002-2003 with the value of J 4 measured 
by GRACE in 2002-2003 (at least within the small measurement uncertainty 
given in the EIGEN-GRACE02S model; this small uncertainty has however 
a negligible influence on our final result). This means that, since we have 
to be consistent with the GRACE observation of J 4 in 2002-2003, the term 
b\ must be: b\ = —2b 2 N, where N is the number of years between the first 
year of application of the J 4 correction and the GRACE observations in 2002- 
2003. Therefore, as we did, one only needs to fit for b 2 and get bi from it, or 
alternatively (as we also did) fit for b 2 (t 2 — 2Nt) in addition to the ait trend. 
So our test fit is just a polynomial fit (with or without the periodic terms). In 
conclusion, in the two sections reported above, the author of (Iorio 2005) is 
explicitly stating that one cannot fit some observational data with a polynomial 
of degree three!! 



2 On the value of j 4 

Long period and secular variations in the zonal harmonics of the spherical 
harmonic expansion describing Earth's gravitational field are generally at- 
tributed to "post-glacial rebound" and ice-sheet mass changes (Trupin, Meier 
and Wahr 1992, Ivins, Sammis and Yoder 1993, Mitrovica and Peltier 1993), 
and reflect changes in Earth's viscous lower mantle (Ivins, Sammis and Yoder 
1993). These changes hold information about the dynamics of Earth and its 
anelastic response to tidal forcing. Estimation of the secular rates for the zonal 
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harmonics therefore has been a high priority research topic of space geodesy 
for many years now, with a considerable number of independent results in 
the literature (Yoder et al. 1983, Rubincam 1984, Eanes et al. 1989, Gegout 
and Cazenave 1991, Nerem et al. 1993, Cheng et al. 1993, Cazenave, Gegout 
and Ferbat 1996, Nerem and Klosko 1996, Eanes and Bettadpur 1996, Cheng, 
Shum and Tapley 1997, Devoti et al. 1999, Cox, Klosko and Chao 2000, Pavlis 
2002, Cheng and Tapley 2004, Cheng and Tapley 2004). The majority of these 
results were obtained from the analysis of long time intervals spanning several 
decades of data obtained from several satellites including the two LAGEOS, 
Starlette, Ajisai, Stella, BE-C and the two ETALONs. Due to the limited 
number of satellite orbit configurations, only a handful of the very lowest de- 
gree zonal terms' rates can be "observed" (Eanes 1995), reaching at best, up 
to degree six. The estimation of the zonal rates is in most cases based on per- 
turbation analysis of very long arcs (several years), but not always necessarily 
so. In particular, the last four studies cited above involve short-arc techniques, 
estimating the average "static" value of the zonal harmonics, and deriving the 
rate from a subsequent analysis of the resulting time series. Furthermore, in 
most cases, the investigators used their static gravitational model of prefer- 
ence, implying that the adopted epoch values for J 4 were different amongst 
studies. These differences are well documented in the literature dealing with 
the static models JGM-3, EGM96, GRIM5-S1, TEG-4, EIGEN-GRACE02S, 
and GGM01S. It is then remarkable that from a plethora of approaches, ana- 
lyzed data sets, and diverse analysis s/w packages used, the resulting estimate 
for the secular rate of these harmonics, e.g. J4, which is of interest here, dif- 
fers by no more than 100% from the commonly accepted value (McCarthy 
and Petit 2004), the one that was also used in the development of EIGEN- 
GRACE02S. 

As it is discussed in (Cheng, Shum and Tapley 1997), there is consider- 
able correlation between the zonal rate estimates (Pj 2 j 4 = — 0-75,Pj 2 j 6 = 
+0.76, pj j = —0.86), but there is no significant correlation between the static 
values of the harmonics and their rates. This independence of errors between 
the static values' estimates and the associated rates is further demonstrated 
in (Cheng and Tapley 2004), where monthly estimates of the J 2 term obtained 
from GRACE data and those obtained entirely independently from SLR data 
analysis are compared. The agreement is remarkable, even though the GRACE 
data analysis has not reached yet its definitive level. The results from (Cheng 
and Tapley 2004) suggest that the current best estimate for J 4 , with all other 
zonal rates up to J 6 simultaneously estimated, is —1.99 x 10~ n /y. This value 
is only 0.58 x 10~ n /y off the currently accepted value, adopted in developing 
EIGEN-GRACE02S, or in terms of a percentage, it is off by only -41%. This 
possible range of error is well within our p error margin (5-10%), since our 
simulation of a 611% error in J4 resulted only in a corresponding 1% error in 
the fitted value of /i. 
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A careful collection of all the published values of J4 in the literature and sub- 
sequent statistical analysis, points to the same conclusion. Since the published 
estimates can be in most cases considered independent as we explained above, 
we can perform a simple statistical analysis of these values and their variance 
in comparison to the quoted error estimates of these values. The population 
is not large by any statistical standards, but that is what is available. There 
are two estimates, +0.2 x 1CT 11 (Nerem and Klosko 1996) and +0.3 x 10 -11 
(Eanes et al. 1989), which are clearly outliers and the reasons are explained in 
(Cheng, Shum and Tapley 1997), so our statistical tests include cases where 
both of these or either one separately are eliminated from the population of 
estimates. The results in terms of a mean J4 and its associated scatter are 
shown in Table [T] 

These tests point to one conclusion very clearly: the uncertainty in our present 
knowledge of J4 is in the neighborhood of ±0.7 x 10 _11 /y and this quantity 
is only 50% of the J4 value that we used in our Lense-Thirring measurement. 
It implies that with a 67% probability, the true value of J4 could be in the 
interval [-1.697 - 0.7 to -1.697 + 0.7] , that is [-2.4, -1.0], or if we increase the 
probability to 95% confidence level, we then have the corresponding interval 
of [-3.1, -0.3], or even at the 99% level, J4 will then lie in the interval [-3.8, 
+0.4], all values in 10~ n /y 

From all these cases, we see that we never reach the extreme valueon: —8.61 x 
10~ n /y, although, even in that case, the final effect on the fitted value of fi 
was only 1%! 

In conclusion, from all the cases studied, there is no reason to expect an error in 
J4 of more than 0.7x 10 _11 /y (at 1-a level of significance, i.e. 67% probability), 
and all statistically acceptable variations in the estimate of J4 indicate that 
there can be no more than 1% error in the fitted value of /x, caused by this error 
in J 4 , as shown in section 1. The evaluation of the J4 error estimate on the 
basis of the various results indicates also a perfect agreement with the quoted 
1 x 10~ n /y error for the —1.41 x 10~ n /y rate for J4, which was adopted in 
developing EIGEN-GRACE02S. 



3 "Imprint" or "Memory" of the Lense-Thirring Effect 

In regard to other conceivable error sources, a possible bias in our measure- 
ment of /1 might be due to some "imprint" or "memory" of the a priori value 
of the Lense-Thirring effect used in the determination of the gravity model 
[see: (Ciufolini et al. 1997)]. The author of (Iorio 2005) is mentioning this po- 
tential source of error in section 3.2.6 "The a priori "memory" effect of the 
Lense Thirring signature on the adopted Earth gravity model on the Earth 
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Table 1 

Statistical significance of the variations between individual J4 estimates. 



Case Description Mean value of J4 Scatter of J4 

(xlO- n /y) (xlO- n /y) 



References 



A. Six estimates 
from Table 5 
(ref. 1) and 

the new estimate 
from (ref. 2) 

B. As in A, 
but without 
the two outliers 
of +0.2 and 
+0.3 x 10- n /y 

C. As in A, 
removing only 
the estimate of 

+0.3 x io- n /y 

D. As in A, 
removing only 
the estimate of 
+0.2 x 10- n /y 

E. Nine estimates 
from Table 7 
of Ref. 1. 



-0.70 



-1.02 



-0.84 



-0.83 



-1.697 



±0.8 (1): (Cheng, Shum and Tapley 1997) 

(2): (Cheng and Tapley 2004) 



±0.6 



±0.74 



± 0.76 



As in (A). 



As in (A). 



As in (A). 



±0.69 (1): (Cheng, Shum and Tapley 1997) 



gravity model used". The Earth gravity field model EIGEN-GRACE02S has 
been obtained using the observations of the GRACE satellites only and by 
using a set of models and parameters to describe the orbital perturbations, 
including some a priori, theoretical, value of the Lense-Thirring effect, indeed 
E1GEN-GRACE02S has been obtained without the Lense-Thirring effect in 
the data analysis. The previous models EIGEN-2S and GGM01S have been 
derived using the observations of the CHAMP and GRACE satellites. Thus, 
the gravitational field model EIGEN-GRACE02S that we used to measure fx 
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may contain some kind of "imprint" of the a priori, theoretical, value of the 
Lense-Thirring effect (i.e., zero in the EIGEN-GRACE02S case). 



Nevertheless, in the present analysis we are only concerned about a conceiv- 
able imprint on the Earth's even zonal harmonics J 2n . The Earth gravitational 
potential is measured by observation of the rate of change of the GRACE 
inter-satellite distance (Reigber et al. 2005), i.e., the gravitational field is 
measured by observation of the relative acceleration of two test particles 
(the equation linking relative acceleration and gravitational field is, in gen- 
eral relativity, the geodesic deviation equation (Ciufolini and Wheeler 1995)). 
However, the effect of the gravitomagnetic field on the acceleration vector 
of a satellite (Ciufolini and Wheeler 1995), decomposed in out-of-plane di- 
rection, _L (orthogonal to the satellite orbital plane); along track, v (along 
the satellite velocity vector); and radial direction, r (along r), is given by 
a = v x H, where H = 2 ^ J ~~ ^ ' x - 1 x ^ is the gravitomagnetic field and v = 

vv. For a polar satellite with a circular orbit we then have a^~ T ~ sine/), 
where is the latitude, a^~ T = and a\~ T = 0. However, the effect of 
the even zonal harmonics on the acceleration vector of a satellite is given by 
a = VU and the terms of the potential U generated by the mass, M, of 
the central body and by the even zonal harmonics, J 2n , are (Kaula 1966): 
U = f[l-J 2 (^) 2 P 2O (sin0) - J 4 (^)4p 4O ( sill 0) _ ...], where J 2 and J 4 are 
the non-normalized even zonal harmonic coefficients of degree 2 (quadrupole) 
and 4, i? e is the Earth's equatorial radius and P 2O (sin0) an d P4o(sin0) are 
the Legendre associated functions. For a polar satellite with a circular orbit 
we then have a± n = 0, a~ 2n ^ and ai 2n ^ 0. Then, for a polar satellite with a 
circular orbit the satellite acceleration generated by the even zonal harmonics 
is just orthogonal to satellite acceleration generated by the gravitomagnetic 
field (i.e. the Lense-Thirring effect). In conclusion, the static even zonal har- 
monics determined by a polar satellite with a circular orbit are independent 
on the a priori value used for the Lense-Thirring effect. 

Indeed, the orbit of the GRACE satellites is nearly circular with eccentricity 
< 0.005 (the same is true for the CHAMP satellite with eccentricity = 0.004) 
and basically polar with inclination = 89° (the same is true for the CHAMP 
satellite with inclination = 87.3°), thus the values of the even zonal harmonics 
determined by the GRACE orbital perturbations are substantially indepen- 
dent on the a priori value of the Lense-Thirring effect. This also applies to 
the CHAMP polar satellite and to the previous CHAMP and GRACE only 
models. The small deviation from a polar orbit of the GRACE satellite, that 
is 1.7 x 10~ 2 rad, gives only rise, at most, to a very small correlation with 
a factor 1.7 x 10~ 2 . For CHAMP this correlation is, at most, ~ 4.7 x 10~ 2 . 
In addition, since our results are independent of the error in the quadrupole 
coefficient J 2 , any bias in J 2 due to a Lense-Thirring "imprint" is canceled in 
our combination (4). However, the Lense-Thirring effect depends on the third 
power of the inverse of the distance from the central body, i.e., (1/r) 3 , and the 
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J2, J±i Je ■■■ effects depend on the powers (1/r) 3 - 5 , (1/r) 5 - 5 , (1/r) 7 - 5 ... of the 
distance; then, since the ratio of the semimajor axes of the GRACE satellites 
to the LAGEOS' satellites is ~ = y^, any conceivable "Lense-Thirring 

imprint" on the spherical harmonics at the GRACE altitude becomes quickly, 
with increasing distance, a negligible effect, especially for higher harmonics of 
degree / > 4. Therefore, any conceivable "Lense-Thirring imprint" is negli- 
gible at the LAGEOS' satellites altitude. The same applies to the CHAMP 
polar satellite. In addition, in (Ciufolini et al. 1997), it was proved with several 
simulations that by far the largest part of this "imprint" effect is absorbed in 
the by far largest coefficient J 2 . 

In conclusion, any error due to a conceivable "imprint" of the Lense-Thirring 
effect on the Earth's even zonal harmonics of the EIGEN-GRACE02S model 
(and of the other CHAMP and GRACE only models) is negligible on the 
combination (4) of the nodes of LAGEOS and LAGEOS II and therefore in 
our measurement of the Lense-Thirring effect. 



4 A brief history of proposed experiments and measurements of 
the Lense-Thirring effect using the nodes of the LAGEOS satel- 
lites 

The measurement reported in Nature (Ciufolini and Pavlis 2004) uses the 
nodes of the two laser ranged satellites LAGEOS and LAGEOS 2 in order to 
cancel the effect of the first even zonal harmonic coefficient, J 2 , of Earth and 
to measure the Lense-Thirring effect. Furthermore, it uses the accurate model 
EIGEN-GRACE02S of the Earth gravity field developed by GFZ of Potsdam 
with the data of the NASA GRACE satellites. 

The key idea of using the nodes of two laser ranged satellites of LAGEOS 
type to measure the Lense-Thirring effect was published for the first time in 
1984-86 (Ciufolini 1984, Ciufolini 1986)). 

The idea to use the nodes of N satellites of LAGEOS type to cancel the effect 
of the first N-l Earth even zonal harmonics and to measure the Lense-Thirring 
effect was published for the first time in 1989 (Ciufolini 1989). The measure- 
ment of the Nature paper is simply the case of N=2. Indeed, in this 1989 
paper on page 3102, fourth line, in order to measure the Lense-Thirring effect 
and to cancel the even zonal harmonics uncertainties, it is proposed (see also 
(Ciufolini and Wheeler 1995), on page 336): "For J 2 , this corresponds from for- 
mula (3.2), to an uncertainty in the nodal precession of 450 milliarcsec/year, 
and similarly for higher Jm coefficients. Therefore the uncertainty in ficiasT 
is more than ten times larger than the Lense-Thirring precession. A solution 
would be to orbit several high-altitude, laser-ranged satellites, similar to LA- 
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GEOS, to measure J 2 , J4, Jq etc, and one satellite to measure Q Lense-Thirring 
" . Thus, the case of the Nature paper is just the one with two satellites. At 
that time the error due to the even zonal harmonics was quite larger due to 
the much less accurate Earth gravity models available at that time and the 
LAGEOS 2 satellite was not yet launched (it was launched in 1992). 

In particular, the idea, to use the two nodes of the satellites LAGEOS and LA- 
GEOS 2 to measure the Lense- Thirring effect, together with the corresponding 
formula also using the perigee of LAGEOS 2, was published for the first time in 
1996 (Ciufolini 1996). See, e.g., the formula (15) on page 1717 of ref. (Ciufolini 
1996): u 5Clj + kiSQjj + k 2 Sujj = /z(31 + 31.5/ci — 57 k 2 ) milliarcsec/ year + 
[other error sources (5Cqo,5Cso, ...,5Ii,5Iii)}" in the slightly different nota- 
tion of the Nature paper this is exactly eq. (1) of (Ciufolini and Pavlis 2004) 

u 01 j 01 r .j 1 • ,1 »t 1 (\Lense-Thirrinq j ALense— Thirring 

where 31 and 31. 5 are written in the JNature paper: iz 7 and iljj 

and SCq , SC 80 , ... are just written £<5J 2n , and k 2 = 0, i.e. there is no use of the 
perigee, but just of the nodes, according to what we suggested in the above 1989 
paper (Ciufolini 1989) and similar following papers. All the relevant LAGEOS 
and LAGEOS 2 nodal rate coefficients of the J 2n , up to J 4 , were explicitly 
given in this 1996 paper. 

In (Peterson 1997) various combinations of the orbital elements of the satellites 
LAGEOS, LAGEOS 2 and LAGEOS 3 were also studied. 

Finally, the combination (1) of ref. (Ciufolini and Pavlis 2004), using the 
nodes only, with the explicit value of the k\ coefficient (that is however an 
absolutely trivial step using the above formula and the coefficients given in 
(Ciufolini 1996), since it corresponds to set k 2 = 0, according to what ex- 
plained in (Ciufolini 1989) and (Ciufolini and Wheeler 1995), and to carry out 
the following highly trivial arithmetics: k\ = —(—9.3/17.2) = 0.54, where the 
coefficients K T 2Q 9* -9.3- 10 11 milliarcsec/yr and K\l 2* 17.2-10 11 milliarcsec/yr 
were explicitly given in (Ciufolini 1996)) was presented and published in the 
SIGRAV 2002 Villa Mondragone (Rome) School proceedings (Ciufolini 2004), 
presented at a plenary talk at the Marcel Grossmann meeeting in Rio de 
Janeiro in July 2003 (to appear in the proceedings (Ciufolini 2003)) and at 
a number of other meetings since 2002. In almost every talk given by one 
of us from 1996 to 2004, it was discussed the possibility of using the nodes 
only of LAGEOS and LAGEOS 2 to just cancel the J 2 effect and measure the 
Lense-Thirring effect. Therefore, the author of (Iorio 2005) seems to have just 
rediscovered in 2003-2004 what already published by one of us many years 
ago. 

The use of the GRACE-derived gravitational models, when available, to mea- 
sure the Lense-Thirring effect with accuracy of a few percent was, since many 
years, a well known possibility to all the researchers in this field and was pre- 
sented by one of us during the SIGRAV 2000 conference (Pavlis 2002) and 
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published in its proceedings, and was published by Ries et al. in the proceed- 
ings of the 1998 William Fairbank conference and of the 2003 13th Int. Laser 
Ranging Workshop (Ries et al. 2003, Ries, Eanes and Tapley 2003), as very 
clearly reported in the Supplementary Discussion of (Ciufolini and Pavlis 2004) 
where it is written: "... (see also ref. [10] by Ries et al.-2003 concluding that, 
in the measurement of the Lense-Thirring effect using the GRACE gravity 
models and the LAGEOS and LAGEOS 2 satellites: "a more current error 
assessment is probably at the few percent level") 

In conclusion, all the claims of the author (Iorio 2005) are simply lacking of 
any rational basis: above is shown how much work was already published on 
this topic before the author of (Iorio 2005) even began to produce any of his 
paper on this topic and to rediscover some earlier results (the papers reported 
above are just a small subset of all the work published by us in the last 20 
years). To avoid the misunderstandings of (Iorio 2005), it would have just 
been a matter of very carefully reading the previously existing literature on 
this subject! 

Let us finally briefly comment on some historical milestones in regard to the 
error analysis of the Nature paper and to the uncertainties in the perturbations 
of the nodes of the LAGEOS satellites. 

The error analysis of our 2004-measurement of the Lense-Thirring effect with 
the LAGEOS and LAGEOS 2 satellites is substantially the same as the error 
analysis of the LAGEOS III/LARES experiment that was carried out in a 
very large number of papers and studies of NASA, ASI and ESA (Tapley, 
Ciufolini et al. 1989-2004). In particular, in 1989, this analysis was performed 
under the supervision of an official NASA committee chaired by Peter Bender, 
with members J. Anderson, J. Armstrong, J. Breakwell, D. Christodoulidis, F. 
Everitt, E. Guinan, R. Hellings, R. King, I. Shapiro, D. Smith and R. Weiss; 
then, in 1994, J. Ries et al., published a revised error budget of the LAGEOS 
III experiment, with a total error of 4 % [this 1994 error analysis of Ries et 
al. is reported with its reference in the Supplementary Discussion of (Ciufolini 
and Pavlis 2004)] and, in 1998, with a total error of 3%-4% [see: (Ries et 
al. 1998) in (Ciufolini 1998)], as reported in the Supplementary Discussion of 
(Ciufolini and Pavlis 2004). 

For differences and similarities between the error analysis of the proposed 
LAGEOS III/LARES experiment and the one of the recent measurement of 
the Lense-Thirring effect with LAGEOS and LAGEOS 2 (Ciufolini and Pavlis 
2004), see the Supplementary Discussion of (Ciufolini and Pavlis 2004) and 
ref. (Ciufolini et al. 2005). Here we just stress that the studies of the LAGEOS 
III/LARES experiment were fundamental milestones in this field, indeed they 
analyzed in details all the main modelling errors in the nodal rates of LAGEOS 
type satellites. 
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5 On the use of the mean anomaly and on the use of Jason to 
measure the Lense-Thirring effect proposed in (Iorio 2004) 



Besides the deep misjudgements and miscalculations of the impact of J 4 and 
of any "Lense-Thirring imprint" in the uncertainty of our measurement of 
frame-dragging, one of the most profound mistakes and misunderstandings of 
(Iorio 2005) is the proposed use of the mean anomaly of a satellite to measure 
the Lense-Thirring effect (in some previous paper by the same author the use 
of the mean anomaly was also explicitly proposed; see e.g. L. Iorio "Some 
comments on the recent results about the measurement of the Lense-Thirring 
effect in the gravitational field of the Earth with the LAGEOS and LAGEOS 
II satellites" arXiv:gr-qc/0411084 v5 19 Apr 2005, however its latest ver- 
sions, after the online publishing of our present paper in New Astronomy on 
the 2nd of May 2005, have been amended of this mistake). In (Iorio 2005) 
(arXiv:gr-qc/0411024 v9 19 Apr 2005, however its latest versions, after the 
online publishing of our present paper in New Astronomy on the 2nd of May 
2005, have been amended of this mistake) is written: 

"The problem of reducing the impact of the mismodelling in the even zonal 
harmonics of the geopotential with the currently existing satellites can be 
coped in the following way. Let us suppose we have at our disposal N (N > 
1) time series of the residuals of those Keplerian orbital elements which are 
affected by the geopotential with secular precessions, i.e. the node, the perigee 
and the mean anomaly: let them be \1/ A , A = LAGEOS, LAGEOS II, etc. Let 
us write explicitly down the expressions of the observed residuals of the rates 
of those elements 8^^ hs in terms of the LenseThirring effect ^f^, of N - 1 
mismodelled classical secular precessions ^f8Ji induced by those even zonal 
harmonics whose impact on the measurement of the gravitomagnetic effect is 
to be reduced and of the remaining mismodelled phenomena A which affect 
the chosen orbital element: 

5^ bs = si>i T + mfsjt + a a , 

A = LAGEOS, LAGEOSII, ... (2) 



The parameter 7 [i is equal to 1 in the General Theory of Relativity and in 
Newtonian mechanics. The coefficients tyf are defined as 

% = 9 -^P- (3) 
oJi 



and have been explicitly worked out for the node and the perigee up to degree 
/ = 20 in Iorio (2002b, 2003a); they depend on some physical parameters of 
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the central mass (GM and the mean equatorial radius R) and on the satellite's 
semimajor axis a, the eccentricity e and the inclination %. We can think about 
Eq. (2) as an algebraic non-homogeneuous linear system of N equations in N 
unknowns which are /i and the N - 1 5Jf. solving it with respect to /i allows 
to obtain a linear combination of orbital residuals which is independent of the 
chosen N - 1 even zonal harmonics." 

This is simply a nonsense paragraph: let us, for example, consider a satellite 
at the LAGEOS altitude, the Lense-Thirring effect on its mean longitude is of 
the order of 2 meters/yr, however the mean longitude change is about 1.8 • 10 11 
meters/yr. Thus, from Kepler's law, the Lense-Thirring effect corresponds to a 
change of the LAGEOS semi-major axis of about 0.09 millimeters. Since, even 
a high altitude satellite such as LAGEOS showed a semimajor axis change of 
the order of 1 millimeter /day, due to atmospheric drag and to the Yarkoski- 
Rubincam effect (because of atmospheric drag, the change of semimajor axis 
and mean motion is obviously much larger for lower altitude satellites), and 
since the present day precision of satellite laser ranging is, even in the case of 
the best SLR stations, of several millimeters, it is a clear nonsense to propose 
a test of the Lense-Thirring effect based on using the mean anomaly of ANY 
satellite, mean anomaly largely affected by non- conservative forces. 

In (Iorio 2005) and in one of Iorio's recent papers, it is proposed to use the 
data collected from the satellite mission JASON, to accurately measure the 
Lense-Thirring effect. Let us highlight some of the major obstacles that one 
would face in trying to implement such a proposal. A more detailed work 
discussing this proposal, the use of the mean anomaly and other highly unfea- 
sible proposals of the author of (Iorio 2005) will be the subjects of a following 
paper. 

First of all, Jason is a "low Earth orbit" satellite, at a 1330 km altitude, and 
experiences by far larger perturbations than the LAGEOS satellites from at- 
mospheric drag, albedo and solar radiation pressure and from all gravitational 
wavelengths. The complicated shape of the satellite makes impossible to pre- 
cisely "forward-model" all the non-conservative forces acting on the satellite 
which generate perturbations that are making JASON a totally unsuitable 
target for making a measurement of such a tiny effect as the Lense-Thirring 
drag. 

JASON is an oceanographic mission which measures radar ranges from the 
satellite to the surface of the ocean with a precision of about 4 cm. To make 
use of these observations for oceanographic research, the requirement on the 
orbital accuracy is stringent but certainly not extraordinary. A couple of cen- 
timeters of radial error are adequate, although recent results indicate that 
we may be closer to 1 cm (Luthcke et al. 2003). The JASON bus is a very 
large, asymmetrical one and it was not designed to support millimeter ranging 



23 



work, since oceanography is still trying to rationalize much larger signals. In 
addition to the large size and poorly defined shape in terms of properties, 
etc., the satellite undergoes nearly continuous orbital adjustments to main- 
tain its ground-track within ± 1 km off the accepted mean ground-track of 
the Topex/Poseidon mission and in this process its ideal orbit as a test particle 
in the Earth gravitational field changes by thrusting and maneuvering; this 
is never easily controlled and it cannot be taken into account in computing 
a perfectly reconstituted orbit. The shape of the satellite bus and the con- 
tinuous motion of parts on the bus (e.g., solar panels) as well as the satellite 
itself (yawing or not, to maintain orientation with respect to the Sun) gener- 
ate further accelerations that are unknown and extremely difficult to account 
for, except by solving for what we call "ad hoc", non-physically meaningful, 
accelerations at certain intervals and that, however, will absorb signals such 
as the L-T effect. To avoid these problems a large international group of re- 
searchers proposed the LARES/WEBERSAT mission (Ciufolini 1998) and, 
in doing so, designed a new structure for the LARES satellite, a very clean 
cannonball satellite, similar to LAGEOS but improved for what regards the 
effect of the non-gravitational perturbations. In any case, one wonders why 
the author of (Iorio 2005) has not proposed the use of which not only have by 
far better shape and orbital stability compared to JASON but they also carry 
ultra precise accelerometers that measure all non-conservative accelerations, 
so that they are in practice "free-falling" particles in vacuum, at least to the 
extent that is covered by the accuracy of these instruments. 

In conclusion, not only all the criticisms in Iorio's paper are completely un- 
founded and misdirected and our total error bubget is fully confirmed to be 
between 5 % and 10 % of the Lense-Thirring effect, but Iorio's paper shows 
a profound lack of understanding and knowledge of the real, practical, data 
analysis of laser ranging observations and of the actual modeling of the Earth 
gravity field. 
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